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CO I In this paper we study the long time behavior of the classical solutions to a 

hydrodynamical system modeling the flow of nematic liquid crystals. This system 
consists of a coupled system of Navier-Stokes equations and kinematic transport 
equations for the molecular orientations. By using a suitable Lojasiewicz-Simon 
type inequality, we prove the convergence of global solutions to single steady states 
as time tends to infinity. Moreover, we provide estimates for the convergence rate. 
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1 Introduction 

We consider the following hydrodynamical system that models the flow of liquid crystal 
materials 

> : 

v t + v ■ Vv-uAv + VP = -AV • [VdQVd+ (Ad - f(d)) <g> d], (1.1) 

V-v = 0, (1.2) 

d t + v Vd -d- Vv = 7(Ad-/(d)), (1.3) 

in Q x (0, oo). Here, Q is a unit square in IR n , (n = 2,3) (the more general case Q = 
il™ =1 (0, Li) with different periods Lj in different directions can be treated in a similar way). 
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v is the velocity field of the flow and d represents the averaged macroscopic/continuum 
molecular orientations in IR n , n = 2, 3. P is a scalar function representing the pressure, 
which includes both the hydrostatic and the induced elastic part from the orientation 
field. The constants u, A and 7 stand for viscosity, the competition between kinetic energy 
and potential energy, and macroscopic elastic relaxation time (Debroah number) for the 
molecular orientation field, respectively. f(d) = ^ (|<i| 2 — l)d with < 77 < 1 may be seen 
as a penalty function to approximate the constraint \d\ = 1, which is due to liquid crystal 
molecules being of similar size [18]. f(d) is the gradient of the scalar valued function 
F(d) = ^drfl 2 — l) 2 . VdOVd denotes the nx n matrix whose (i, j)-th entry is given by 
Vid ■ Vjd, for 1 < i, j < n. ® is the usual Kronecker multiplication, e.g., (a ® b)ij = afij 
for a, be K™. 

The above system is a simplified version of the Ericksen-Leslie model for the hydro- 
dynamics of nematic liquid crystals (cf. [2,3,8,15,16]). Generally speaking, the system is 
a macroscopic continuum description of the time evolutions of these materials influenced 
by both the flow field v(x,t), and the microscopic orientational configuration d(x,t), 
which can be derived from the coarse graining of the directions of rod-like liquid crys- 
tal molecules. Equation (11.11) is the conservation of linear momentum (the force balance 
equation). It combines a usual equation describing the flow of an isotropic fluid with 
an extra nonlinear coupling term, which is anisotropic. This extra term is the induced 
elastic stress from the elastic energy through the transport, represented by the third equa- 
tion. Equation (11.21) represents incompressibility of the fluid. Equation (11.31) is associated 
with conservation of the angular momentum. The left hand side of (11.31) stands for the 
kinematic transport by the flow field, while the right hand side represents the internal 
relaxation due to the elastic energy. The continuum theory of liquid crystals due to Er- 
icksen and Leslie was developed during the period of 1958 through 1968. Since then there 
has been a remarkable research in liquid crystals, both theoretically and experimentally. 

In the context of hydrodynamics, the basic variable is the flow map (particle trajectory) 
x(X, t). X is the original labeling (the Lagrangian coordinate) of the particle, which is 
also referred as the material coordinate, x is the current (Eulerian) coordinate, and is also 
called the reference coordinate. For a given velocity field v(x, t), the flow map is defined 
by the ODE : 

x t = v(x(X,t),t), x(X,0) = X. 

To incorporate the elastic properties of the material, we need to introduce the deformation 
tensor 

HX,t) = —(X,t). 
This quantity is defined in the Lagrangian material coordinate and it satisfies 

dF(X, t) _ dv(x(X,t),t) 
dt ~ dX ' 
In Eulerian coordinates, we define T{x, t) = J-{X, t). We shall use the notation Tij = 
^jrh. By using the chain rule, the above equation can be transformed into the following 
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transport equation for T (cf. [6]): 

F t + (vV)r = VvF. 

Without ambiguity, in the following text, we will not distinguish these two notations T 
and T . If the liquid crystal has a rod-like shape, then transport of the direction of d can 
be expressed as 

d(x(X,t),t) =Fd Q {X), 

where do(X) is the initial condition. This equation demonstrates the stretching of the 
director besides the transport along the trajectory. By taking full time derivative on both 
sides, we have (cf. [26]) 

^-d(x(X, t),t) = Fd {X) = VvFdo = Vvd = (d ■ V)v. 
Hence, the total transport of the orientation vector d becomes 

dt + v • Vd — d ■ Vf , 

which represents the covariant parallel transport with no-slip boundary condition between 
the rod-like particle and the fluid (cf. [10,14]). In general, for a molecule of ellipsoidal 
shape with a finite aspect ratio, the transport of the main axis direction is represented by 

d(x{X,t),t) = Edo(X), 

where E is a linear combination of T and T~ T and satisfies the transport equation: 

E t + (v V)E = (aVv + (1 - a)(-V T v))E. 

The spherical, rod-like and disc-like liquid crystal molecules correspond to a = ~, 1 and 
0, respectively (cf. [1,2,19,26]). 

In this paper, we shall focus on the rod-like nematic liquid crystal molecules. Analysis 
for the case of general ellipsoidal molecules may be done in the future work. The system 
(jl.lD - fll.3p can be derived from an energetic variational point of view. Consider the 
following action functional in terms of flow map x(X,t): 



A 



o Jn 



i|x t | 2 - \ l -\F~ T V x Td {X)\ 2 - \F(Fd (X)) 



JdXdt, 



where f2 is the region at the initial time and J = detj^7 is the Jacobian. It has been 
shown that by using the least action principle (Hamilton's principle), 4^ = 0, we can 
recover the system fll . IP — fll .3H . We refer to [26] for details. It has been pointed that 
such a derivation using the variation with respect to domain, i.e., least action principle, 
is equivalent to the principle of virtual work. 
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In this paper, we will consider system fll .lf> — fll .3D subject to the periodic boundary 
conditions 

v (x + ej) = v(x), d(x + ej) = d(x), for x G dQ, (1.4) 
and initial conditions 

v\t=o = vo(x) with V ■ t>o = 0, d\ t=Q = d (x), for x e Q. (1.5) 

Here, {i = 1, n) are the canonical basis of W 1 . The global existence of weak/classical 
solutions to the system (11.1I) - (11.5I) for n = 2 or n = 3 with large viscosity assumption has 
been proven in [26]. It was pointed out that by supposing periodic boundary conditions, 
we can get rid of the boundary terms when integrating by parts, which is necessary in the 
derivation of higher-order energy inequalities. Without this assumption, some boundary 
terms would remain persistent and undermine the whole estimates. The corresponding 
initial boundary value problems are still open. 
The main results of this paper are as follows 

Theorem 1.1. When n = 2, for any initial data (t>o, do) G V x H^{Q), the unique global 
classical solution to problem fll . IB — fll .5D has the following property: 

lim (\\v(t)\\ m + \\d(t)-d 00 \\ H 2)=U } (1.6) 

t— >+oo 

where doo is a solution to the following nonlinear elliptic boundary value problem: 

- Ad OQ + f(d oc ) = 0, xeQ, d^x + ei) = d^x), x e dQ. (1.7) 

Moreover, there exists a positive constant C depending on vq, d , Q, such that 

|Kt)||tfi + |Kt)-4o||if2 <C(l + t) _ ^T, Vt>0. (1.8) 

Here, 6 e (0, 1/2) is the same constant as in the Lojasiewicz-Simon inequality (see Lemma 
[PI below). 

In three dimensional case, we have the following results: 

Theorem 1.2. When n = 3, for any initial data (vo,do) G V x Hp(Q), under the large 
viscosity assumption v > Vq(X, 7, Vq, do) (cf. (14.24ft ). the problem (11.11) - (11.51) admits a 
unique global classical solution enjoying the same properties as in Theorem \l.l[ 

Theorem 1.3. When n = 3, let d* G H^[Q) be an absolute minimizer of the functional 

E(d) = h\Vd\\ 2 + [ F{d)dx (1.9) 
2 Jn 

in the sense thatE(d*) < E(d) for alld G H^{Q). There exists a constant a G (0, 1], which 
may depend on A, 7, /, Q and d*, such that for any initial data (vo,d ) G V x H^{Q) 
that satisfy ||f ||^i + ||c?o — ^*||h" 2 < °~ , problem fll.ip - fll.5D admits a unique global classical 
solution enjoying the same properties as in Theorem \l.l\ 
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Remark 1.1. Theorem \1.3\ implies that if the initial velocity Vq is small and the initial 
molecule orientation oIq is sufficiently close to an absolute minimizer of the functional E, 
then problem (11.11) - (11.51) admits a unique global solution and the solution will converge to 
a certain equilibrium, which is not necessarily the absolute minimizer itself. 

Theorems 1.1-1.3 imply the uniqueness of asymptotic limit of system (II. IB - Ill. 5ft when- 
ever it admits a global bounded solution. The problem on uniqueness of asymptotic limit 
for nonlinear evolution equations, namely whether the global solution will converge to an 
equilibrium as time tends to infinity, has attracted a lot of interests among mathemati- 
cians for a long time. If the spacial dimension n > 2, it is known that the structure of 
the set of equilibria can be nontrivial and may form a continuum for certain physically 
reasonable nonlinearities. In particular, in our present case, under the periodic boundary 
conditions, one may expect that the dimension of the set of stationary solutions is at least 
n. This is because a shift in each variable should give another steady state (cf. [23]), e.g., 
if d* is a steady state, so is d*(- + re^), 1 < i < n, r G R + . In this case it is highly 
nontrivial to decide whether a given trajectory will converge to a single steady state. In 
1983, L. Simon [25] made a breakthrough by proving that for a semilinear parabolic equa- 
tion with a nonlinear term f(x,u) being analytic in the unknown function u, its bounded 
global solution would converge to an equilibrium solution as t — > <x>. Simon's idea re- 
lies on a generalization of the Lojasiewicz inequality (cf. [21,22]) for analytic functions 
defined in the finite dimensional space M. m . Since then, his original approach has been 
simplified and applied to prove convergence results for various evolution equations (see 
e.g., [4,5,7,9,11,30] and the references therein). In order to apply the Lojasiewicz-Simon 
approach to our problem (ll.ip -f Tl~5l) . we need to introduce a suitable Lojasiewicz-Simon 
type inequality for vector functions with periodic boundary condition (cf. Lemma [2.3ft . 
The corresponding result for small molecule system (cf. [18]) was shown in [20] under var- 
ious boundary conditions (Dirichlet b.c. /free-slip b.c.) and also for flows with changing 
fluid density. 

Once we prove the convergence to an equilibrium, a natural subject for further study 
is the convergence rate. It is well known that an estimate in certain (lower-order) norm 
can usually be obtained directly from the Lojasiewicz-Simon approach (see, e.g., [7,31]). 
One can then in a straightforward way, obtains estimates in higher-order norms by using 
interpolation inequalities (cf. [7]), and consequently, the decay exponent deteriorates. We 
shall show that by using suitable energy estimates and constructing proper differential 
inequalities, it is possible to obtain the same estimates on the convergence rate in both 
higher and lower order norms. Our approach improves the previous results in the literature 
and can be applied to many other problems (cf. [5,29,30]). 

The remaining part of this paper is organized as follows. In Section 2, we introduce 
the functional setting, some preliminary results and some technical lemmas. Section 3 
is devoted to the two dimensional case. We prove the convergence of global solutions 
to single steady states as time tends to infinity and obtain an estimate on convergence 
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rate. In Section 4, we consider the three dimensional case. The same convergence result 
is proved for two subcases, in which the global existence of classical solution can be 
obtained. In the finally Section 5, we briefly discuss the results for liquid crystal flows 
with non-vanishing average velocity. 

2 Preliminaries 

Recall the well-established functional setting for periodic problems (cf. [18,27]): 

H'p{Q) = {v G H m (R n ) I v(x + et) = v(x)}, 
H™(Q) = H;\Q) n : J v(x)dx = | , 

H = {v G Lj(Q), V • v = 0}, where L 2 p (Q) = H° p {Q), 

V = {veH^Q), V-v = 0}, 

V = the dual of V. 

For simplicity, we denote the inner product on L p (Q) as well as H by (-, •) and the 
associated norm by || ■ ||. We shall denote by C the genetic constants depending on 
A, 7, v, Q, f and the initial data. Special dependence will be pointed out explicitly in the 
text if necessary. Since the parameters A and 7 do not play important roles in the proof, 
we set A = 7 = 1 for the sake of simplicity. The following embedding inequalities will be 
frequently used in the subsequent proofs: 

Lemma 2.1. (cf. [27]) If n = 2, we have 

1 1 

IMU»(Q) < c IMMMIJr2(Q)> vmg h%(q), 

Ifn = 3, then 

1 — o 

IMIl-(q) < c|MH|w||^ 2(Q) , V u e H p (Q). 

Here, we note that \\u\\h 2 {q) can be estimated by ||Au|| and \\u\\ in sprit of the elliptic 
estimate ( 12.24B . 

The global existence of weak/classical solutions to system p.ip -f fl~5l ) has been proven 
in [26, Theorem 1.1]. More precisely, we have 

Proposition 2.1. Assume that (va,do) G V x H p {Q). Then, if either n = 2 or n = 3 
with the large viscosity assumption v > C(X, 7, Vq, do) (see (14.241) ). problem (ll.ip - (ll.5l) 
admits a global solution such that 

v e L°°(0,oo;V), d G L°°(0, 00; H 2 ), (2.1) 

and for any T > 0, 

v G L 2 (0, T; H 2 ), d G L 2 (0, T; if 3 ). (2.2) 

Moreover, v,d G C°°(Q) for all t G (0,T), (v,d) and all of their spatial derivatives are 
absolutely continuous in time. 
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The proof for Proposition 12. II relies on a modified Galerkin method introduced in [18]. 
After generating a sequence of approximate solutions, one can use the Ladyzhenskaya 
method (cf. [13,28]) to get a high-order energy estimate, which enables us to pass to the 
limit. Furthermore, a weak solution together with high-order derivative estimates implies 
a strong solution, i.e. v G L 2 (0,T; H 2 (Q)) and d G L 2 (0, T; H 3 (Q)). Finally, a bootstrap 
argument based on Serrin's result [24] (cf. also [12]) and Sobolev embedding theorems 
leads to the existence of classical solutions. Comparing with the small molecule system 
(cf. [18]), we now have different kinematic transport and accordingly one more stress term 
(Ad — f(d)) Cg) d in the elastic stress in (II. ID and one more transport term d ■ Vt> in fll .3D . 
These bring extra technical difficulties to prove the existence result. For instance, d ■ Vv 
stands for the parallel transport, which includes both rotation and stretching effect of the 
director d. The stretching effect leads to the loss of maximum principle for the equation 
for d. On the other hand, the extra stress term (Ad — f(d)) ® d cannot be suitably defined 
in the weak formulation as in [18], and thus the requirement that d G L°°(0, T; L°°(Q)) 
must be imposed so that the problem is well-posed. For detailed discussions, we refer 
to [26]. 

The Lyapunov functional for problem (ll.ll) - (ll.5l) is 

S(t) = l -\\v(t)\\ 2 + ^||W(i)|| 2 + A ! F(d(t))dx. (2.3) 
Our system has the following basic energy law 

j t £(t) = -v\\Vv(t)\\ 2 - A 7 ||Ad(f) - f(d(t))\\ 2 , t > 0. (2.4) 

It is worth pointing out that, for different molecule shapes (and as a result, different 
kinematic transports), the system possesses the same energy dissipative law (12.41) (cf. [19]). 

First, we shall show a continuous dependence result on initial data, from which the 
uniqueness of classical solutions to problem (II. IB - Ill. 5D follows. 

Lemma 2.2. Suppose that (vi,di) (i = 1,2) are global solutions to problem (11.11) - (11.51 ) 

corresponding to initial data (voi,d 0i ) G V x H 2 (Q) (i = 1,2). Moreover, we assume that 
for any T > 0, the following estimate holds 

IK(*)||iji + R(t)||tf 2 <M, VtG[0,T]. (2.5) 

Then for any t G [0, T], we have 

\\{vi - v 2 )(t)\\ 2 + Ufa - d 2 )(t)\\m + /Vl|V(«i - v 2 )(t)\\ 2 + ||A(d! - d 2 )(r)\\ 2 )dr 
< 2e ct (\\v Ql -v Q2 \\ 2 + \\d m -d Q2 \\ 2 m ), (2.6) 

where C is a constant depending on M but not on t. 
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Proof. Denote 



v = Vi — v 2 , d = d\ — d 2 . (2.7) 

Since (uj,dj) are solutions to problem fll.lH — fll.5D . we have 

v u + vi ■ Vvi - vA Vl + VPi = -V-[Vdi0Vdi + (AtZi-/(di))®di], (2.8) 

V-vi = 0, (2.9) 

dit + «i • Vdi - di • Vui = Ad x -/(di), (2.10) 

and 

u 2t + u 2 • Vu 2 - z/Au 2 + VP 2 = -V- [Vd 2 © Vd 2 + (Ad 2 -/(d 2 ))®d 2 ], (2.11) 

V-u 2 = 0, (2.12) 

d 2t + v 2 -Vd 2 -d 2 -Vv 2 = Adz- f ((h) . (2.13) 

Multiplying Ui — v 2 with the subtraction of (12.111) from (12.81) . and multiplying (di — d 2 ) — 
(Adx — Ad 2 ) with the subtraction of (12.131) from (12.101) . we add these two resultants 
together. By integration by parts, we infer from the periodic boundary conditions that 

\j t (IN 2 + W + livdf ) + Hiv^ll 2 + ||vd|| 2 + ||AJ|| 2 

= -(u 2 ■ Vu, u) - (u • Vu x , u) - (Ad 2 ■ Vd, v) + (Ad 2 ® d, Vu) 

-((/(di) - /(d 2 )) di, Vu) - (f(d 2 ) d, Vu) + (M) - /(d 2 ), Ad) 

+(i; 2 ■ Vd, Ad) - (d ■ Vu 2 , Ad) - (/(di) - /(d 2 ), d) - (u ■ Vd 1; d) 

-(u 2 ■ Vd, d) + (di • Vu, d) + (d • Vu 2 , d). (2.14) 

From assumption (12 .5p . ||u||#i and ||d||#2 are uniformly bounded in [0, T]. Hence, by 
using the Sobolev embedding theorems, we can estimate the righthand side term by term 
(the calculation presented here is for the three dimensional case and it is also valid for 
two dimensional case). 

|(u 2 - Vu,u)| + |(u- Vu x ,u)| 

< ||^2||l 6 I|V^II ||u||z,3 + ||u|||4||Vui|| 

< C||Vu||(||Vu||i||u||i + ||u||) + C(||Vu||t||w||i + \\v\\) 2 

< e\\Vv\\ 2 + C\\v\\ 2 . (2.15) 



| ( Ad 2 • Vd, u) | + | ( Ad 2 <g> d, Vu) | 

< ||Ad2||||Vdl|rs||t?|U6 + ||Ad 2 ||||d|| L oo||Vu|| 

< C(||Adl|l||Vdl|l + ||Vd||)(||Vu|| + ||u||) + C(||Adl|i + ||d||)||Vu|| 

< e(\\Ad\\ 2 +\\V^ 2 )+C(\\d\\ 2 Hl + \\v\\ 2 ). (2.16) 
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\((f(d x ) - f(d 2 )) <g> d u Vv)\ + \(f(d 2 ) ® d,Vv)\ + |(/(dO - /(d 2 ), Ad)| 
+\(f(d 1 )-f(d 2 ),d)\ 

< (HM) - + ll/WllL-lldlDUVtJii + - f(d 2 )\\(\\Ad\\ + ||J||) 

< C(||/'(OIU~ + l)||dI|||VT?|| +C||/'(OIU-||d||(||Ad|| + ||d||) 

< e(\\Vv\\ 2 + \\Ad\\ 2 ) + C\\d\\ 2 , (2.17) 
where £ = ad\ + (1 — a)d 2 with a G [0, 1]. 

\(v 2 - Vd,Ad)\ + | (J- Vv 2 ,Ad)\ 

< ||w 2 || L e||vJ|| L 3||Aj|| + ||Vw 2 ||||J|| L oc||AJ|| 

< C(||Ad||i||Vd||& + ||Vd||)||Aj|| +C(||Ad||f + H|)||AJ|| 

< e||Ad|| 2 + C||d||^i. (2.18) 



\(v-Vd u d)\ + \(v 2 -Vd, d)\ + \{dx- Vv,d)\ + \(d- Vv 2 , d)\ 

< ||Vdi|U8||i;||||dl| Le + ||^|| i6 ||VcJ||||J|| L:i + HlU-UV^IIIIdll + ||V^ 2 ||||J||1,4 

< e\\Vv\\ 2 + C{\\v\\ 2 +\\df m ). (2.19) 

Choosing e small enough in the above estimates, we infer from (12.141) that 

j t {M\ 2 + \\d\\ 2 m ) + v\\Vv\\ 2 + ||Ad|| 2 < C(|M| 2 + || d\\ 2 Hl ), (2.20) 

where C is a constant depending on ||fj||_f/i and ||g?;||#2 but not on t. By Gronwall's 
inequality, we can see that for any t G [0, T] , 

ll^)H 2 + \\d(t)\\m + f ("\\Vv(t)\\ 2 + l|AJ(r)|| 2 )rfr < 2e c \\\v(0) || 2 + || d(0)\\ 2 m ). (2-21) 
Jo 

The proof is complete. □ 

Remark 2.1. Since the global classical solution (v,d) to problem fll . lh f jT75T ) obtained in 
Proposition ^. 1\ is uniformly bounded in H 1 x if 2 ( cf. [26], see also Lemma lSTB and Lemma 
4-l\ below), it immediately follows from Lemma lKE that this global solution is unique. 



Next, we look at the following elliptic problem with periodic boundary condition 

-Ad + /(d) = 0, xeQ, 



Define 



(2.22) 

d(x + e$) = d(x), x G dQ. 

E{d):=\\\Vd\\ 2 + f F{d)dx. (2.23) 
2 Jq 

It is not difficult to see that the solution to (12.221) is a critical point of E(d) and conversely 
the critical point of E(d) is a solution to (12.221) . Moreover, the solution to (12.221) is 
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smooth. We recall the interior elliptic estimate, which states that for any U\ CC U 2 there 
is a constant C > depending only on U\ and U 2 such that || rfH^r^rc/i) < C(II^^IIl 2 (c/ 2 ) + 
||c?||l2( [ / 2 )). In our case, we can choose Q' to be the union of Q and its neighborhood 
copies. Then we have 

\\d\\ H2{Q) < C(\\Ad\\ LHQ/) + \\d\\ L2m ) = 9C(||Ad|| L2(Q) + \\d\\ LHQ) ). (2.24) 

In order to apply the Lojasiewicz-Simon approach to prove the convergence to equi- 
librium, we shall introduce a suitable Lojasiewicz-Simon type inequality that is related 
to our problem. In particular, we have 

Lemma 2.3. [Lojasiewicz-Simon inequality] Let ifj be a critical point of E(d). Then 
there exist constants 9 £ (0, |) and (3 > depending on ip such that for any d £ H^(Q) 
satisfying \\d — ip\\m(Q) < P, it holds 

|| - Ad + f(d)\\ mm > \E{d) - Em 1 - 9 , (2.25) 

where (Hp(Q))' is the dual space of Hp(Q). 

Remark 2.2. Lemma Wlh can be viewed as an extended version of Simon's result [25] for 
scalar functions using the L 2 -norm. For the proof of this result, we refer to [9, Chapter 
2, Theorem 5.2]. 

The following lemma turns out to useful in the study of long-time behavior of solutions 
to evolution equations. 

Lemma 2.4. [31, Lemma 6.2.1] Let T be given with < T < +oo. Suppose that y(t) and 
h(t) are nonnegative continuous functions defined on [0, T\, which satisfy the following 
conditions: 

^ <c lV 2 + c 2 + h{t), j y{t)dt<c 3 , j h{t)dt<c 4 , 

where Ci (i = 1,2,3,4) are given nonnegative constants. Then for any r £ (0,T), the 
following estimate holds: 

y(t + r) < (^ + c 2 r + c 4 ) e clC3 , Vt £ [0,T - r). 

Furthermore, ifT = +oo, then 

lim y(t) = 0. 

3 Long-time Behavior in Two Dimensional Case 

In this section, we prove the convergence of global solutions to single steady states as time 
tends to infinity in the two dimensional case. In 2-D case, an important property for the 
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global solution to problem fll.lH — fll.5H is the following high-order energy inequality, which 
played a crucial role in the proof of global existence result in [26]. 
Denote 

A(t) = \\Vv(t)\\ 2 + \\\Ad(t)-f(d(t))\\ 2 . (3.1) 

We recall that it has been assumed that A = 7 = 1. Besides, since viscosity v does 
not play a crucial role in the 2-D case, we also set v — 1 in this section for the sake of 
simplicity. Then we have (cf. [26, (45)]) 

Lemma 3.1. In two dimensional case, the following inequality holds for the classical 
solution (v,d) to problem ( ll.ll) -( TIT5"l ): 

jA{t) + \{\\^v{tW + ||V(Ad(*) - /K*)))ll 2 ) < C{A\t) + 1), V t > 0, (3.2) 
where C is a constant depending on f,Q, ||fo||, II^oIIh 1 ^)- 

3.1 Convergence to Equilibrium 

Based on the high-order energy inequality (13 .2D . we are able to show the decay property 
of the velocity field v. 

Lemma 3.2. For any t > 0, the following uniform estimate holds 

\Ht)\\ H i + \\d(t)\\ H i<C, (3.3) 
where C is a constant depending on f,Q, \\voWh 1 , \\do\\H 2 - Furthermore, 

lim (\\v(t)\\ m + || - Ad(t) + f(d(t))\\) = 0. (3.4) 

t— >+oo 

Proof. It follows from the basic energy law (I2.3p that 

£(t) + [ A(r)dr < £ (0) < 00, V t > 0. (3.5) 
Jo 

(13 .5p implies the uniform estimate 

IK*)ll + IK*)llffi<C, Vt>0. (3.6) 

Concerning the uniform bound (13.31) . we take r = 1 in Lemma T2.4I to get 

\\Vv(t) || + || - Ad(t) + f(d(t))\\ <C, V t > 1, (3.7) 

where C does not depend on t. On the other hand, for any t e [0, 1], it follows from (13.21) 
and the fact J* A(t)dt < C that 

sup A{t) < eti A{t)dt A{0) + C <C. (3.8) 

0<t<l 
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Besides, from the continuous embedding H 1 LP (1 < p < oo) and (13.61) we have 

||Ad|| < || -Ad + /(d)|| + \\f(d)\\ < || - Ad + /(d)|| +C(1+ ||d||| 6 ) < C. (3.9) 

Therefore, (I3.3P follows from (I3.7I) - (I3.9I) . Furthermore, (13 .51) together with Lemma [37T1 
and Lemma E31 yields that 

lim (||V V (t)|| + ||Ad(t)-/(d(t))||) = 0. (3.10) 

By the Poincare inequality for v G V, we conclude (13.41) . The proof is complete. □ 
Let S be the set 

S = {(0, «) | — Au + f{u) = 0, in Q, u(x + e$) = m(x) on <9Q}- 

The cj-limit set of {v ,d ) eV x Hp(Q) C x H^(Q) is defined as follows: 

uj(vo,do) = {(v 00 (x),d 00 (x)) | there exists {t n } /* oo such that 

(v(t n ),d(t n )) -> (^00,^00) in L 2 (Q) x H X {Q), as i n -> +oo}. 

We infer from Lemma T3.2I that 

Proposition 3.1. o;(u ,do) is a nonempty bounded subset in Hp(Q) x H^{Q), which is 
compact in Lt\Q) x HHQ). Besides, all asymptotic limiting points ('Uoojdoo) of problem 
(ll.ip - (ll.5p belong to S. In other words, u(vo,d ) C 5. 

In what follows, we prove the convergence for director field d. For any initial data 
(vo, do) € V x Hp(Q), it follows from Lemma [3721 that ||d||#2 is uniformly bounded. 
Proposition 13.11 implies that there is an increasing unbounded sequence {in}ngN and a 
function d^ such that 

lim \\d(t n ) -doolie = 0. (3.11) 
Moreover, d^ satisfies the equation 

- Adoo + /(doo) = 0, a; e O, d^x + e*) = d^x) on 9Q. (3-12) 

We prove the convergence result following a simple argument first introduced in [11] 
in which the key observation is that after a certain time to, d(t) will fall into a certain 
small neighborhood of d^ and stay there forever. 

From the basic energy law (12.31) . we can see that £{t) is decreasing on [0, oo), and it 
has a finite limit as time goes to infinity because it is nonegative. Therefore, it follows 
from (E3D and (I3TTD that 

lim £{t n ) = E{d OB ). (3.13) 

On the other hand, we can infer from (I2.3P that £{t) > E(d oa ), for all t > 0, and the 
equal sign holds if and only if, for alH > 0, v = and d solves problem (13.121) . 
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We now consider all possibilities. 
Case 1. If there is a t > such that 8 (to) = E(d 00 ), then for all t > t , we deduce from 
(Q that 

||Vv||=0, || -Ad + f(d) || = 0. (3.14) 

It follows from ( 11.31 ). ( 13.141) and the Sobolev embedding theorem that for t > to, \\d t \\ = 0. 
Namely, d is independent of time for all t >t . Due to ( 13.111) . we conclude that d(t) = d^ 
for t > t - 

Case 2. For all t > 0, we suppose that £(t) > E(d 00 ). First we assume that the following 
claim holds true. 

Proposition 3.2. There is at > snc/i that for allt > £ , ||d(t) — rfooll// 1 < Namely, 
for all t > t , d(t) satisfies the condition in Lemma \2.3[ 



In this case, it follows from Lemma [2731 that 

\E(d) - ^(cUI 1 - 9 < || - Ad+ f(d)\\ m y < || - Ad+ /(d) II, V t > t . (3.15) 
The fact 9 G (0, ~) implies that < 1 - 9 < 1, 2(1 - 9) > 1. Consequently, 

|| u || 2 ( 1 - e ) = || u || 2 ( 1 - e )- 1 || v ||<C7||i;||. 
Then we infer from the basic inequality 

(a + 6) 1 - < a^ + fe 1 - , Va,6>0 

that 



(^(t)-^)) 1 - 9 < (^\\ v f + \E(d)-E(d 00 )\ S j 



< QiHi 2 +n-Ad+/(d)iii^y e 

< e ||^|| 2(1 - e) + ||-Ad + /(d)|| 

< C|H| + ||-Ad + /(d)||. (3.16) 



Therefore, a direct calculation yields that 

- j t (8(t) - E(d ao )) 6 = -9(8(t)-E(d 00 )) s ^j t 8(t) 

C9(\\Vv\\ + \\-Ad + f(d)\\f 



C\\v\\ + || - Ad + /(d) || 
> C 1 (\\Vv\\ + \\-Ad + f(d)\\) 1 Vt>t , (3.17) 
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where C\ is a constant depending on v , d , Q and 8. 
Integrating from t to t, we get 

(S(t) - Eid^Y + d A||V^(r)|| + || - Ad(r) + f(d(r))\\)dr 

J to 

< (S(t ) - Eid^)) 9 < oo, Vt>t - (3.18) 
Since £(£) — E(doo) > 0, we conclude that 

(||V^(r)|| + || - Ad(r) + f(d(r))\\)dr < oo. (3.19) 



On the other hand, it follows from equation (11.31) . ( I3.3p and Sobolev embedding theorems 
that 

\\d t \\ < \\v-Vd\\ + \\d-Vv\\ + \\-Ad + f(d)\\ 

< \\v\\ L 4Vd\\ L * + ||d|U-||Vt;|| + || - Ad+ /(d) || 

< C7||Vv|| + ||-Ad + /(d)||. (3.20) 



Hence, 



|dt(r)||dr < +oo, (3.21) 



which easily implies that as t — > +oo, d(£) converges in L 2 (Q). This and <RTTT]> indicate 
that 

lim ||d(t) - dooll = 0. (3.22) 

Since d(t) is uniformly bounded in H 2 (Q) (cf. (13.31) ). by standard interpolation inequality 
we have 

lim || d{t) -d^Wffi = 0. (3.23) 

t— >+oo 

On the other hand, the uniform bound of d in H 2 (Q) implies the weak convergence 

d(t) - d^, in H 2 {Q). (3.24) 

However, the decay property of the quantity A(t) (cf. Lemma I3T21 could tell us more. 
Namely, we could get strong convergence of d in H 2 . To see this, we keep in mind that 

HAd-Adooll < \\M-M 00 -f{d)+f{d 00 )\\ + \\f{d)-f{d 00 )\\ 

< ||Ad-/(d)|| + ||f(0|M|d-doo||L4 

< ||Ad-/(d)|| + C||d-d 0O || H a. (3.25) 
The above estimate together with ( 13 A\\ and (13.231) yields 

lim \\d(t) - d^Wjp = 0. (3.26) 
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To finish the proof, it remains to show that Proposition 13.21 always holds true for the 
global solution d(t) to system fll.lH — fll.5D . Define 

t n = sup{ t > t n \ \\d(s) - d^Wni <P, V s G [t n ,t]}. (3.27) 

It follows from (13. 1 1H that for any e G (0, /?), there exists an integer N such that when 

n > N, 

\\d(t n ) - d^Wni < e, (3.28) 
^r{S{t n ) - Eid^)) 6 < e. (3.29) 

On the other hand, we know that the orbit of the classical solution d is continuous in H 1 . 
It follows from (Q that d G L°°(0, +oo; H 2 ). As a consequence, d G L 2 (t,t + 1; H 2 ) for 
any t > 0. The basic energy law and (13.201) imply d t G L 2 (t,t + 1;L 2 ). Thus, for any 
t > 0, it holds d G C([t,t + 1]; i^ 1 ). The continuity of the orbit of d in if 1 and (13T28D yield 
that 

i n > tn, for all n> N. (3.30) 

Then there are two possibilities: 

(i) . If there exists n > N such that i no = +oo, then from the previous discussions in 
Case 1 and Case 2, the theorem is proved. 

(ii) Otherwise, for all n > N, we have t n < i n < +oo, and for all t G [t n ,i n ], 
E(d oa ) < £{t). Then from (13.181) with t being replaced by t n , and t being replaced by t n , 
we obtain from (13.291) that 

t (\\Vv(t)\\ + || - Ad(r) + /(d(r))||)dr < 5. (3.31) 
Thus, it follows that (cf. (13T20D ) 

/•in 

||d(t n ) -dooll < ||d(t n ) - dooll + / \\dt(r)\\dr 



< MQ-d^W+C (||Vu(r)|| + ||-Ad(r) + /(d(r))||)dr 

< Ce, (3.32) 

which implies that lim„^ +00 ||d(t n ) — doo|| = 0. Since d(t) is relatively compact in H 1 , 
there exists a subsequence of {d(i n )}, still denoted by {d{t n )} converging to d^ in H 1 , 
i.e., when n is sufficiently large, 

||d(t n ) - doollni < A 

which contradicts the definition of i n that ||d(t„) — doollff 1 — P- 

Summing up, we have considered all the possible cases and the conclusion (11.6P is 
proved. □ 
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3.2 Convergence Rate 

In this part, we shall prove the estimate for convergence rate (11 .81) . This can be achieved 
by several steps. 

Step 1. As has been shown in the literature (cf. [7,31]), an estimate on the convergence 
rate in certain lower-order norm could be obtained directly from the Lojasiewicz-Simon 
approach. From Lemma [2.31 and (13.171) . we have 

j t (£(t) - Eidoo)) + Ci(£(*) - £(cU) 2(1 ~ e) < 0, V t > t , (3.33) 
which implies 

£{t)-E{d 00 )<C{l + ty^, Vt>t . (3.34) 
Integrating A3. 1 7[) on (t, oo), where t > t , it follows from (13.201) that 

/oo poo 
\\d t (r)\\dT< J (C\\Vv(r)\\ + || - Arf(r) + f(d(r))\\)dr < C{\ + t)~^ . (3.35) 

By adjusting the constant C properly, we obtain 

\\d{t) -dooW < C(l + t)~^, t>0. (3.36) 

Step 2. In Step 1, we only obtain the convergence rate of d (in L 2 ). Unlike for the 
temperature variable in some phase-field systems (cf. [30] and references cited therein), 
although we have got some decay information for the velocity field v such that 

Vv(r)\\dr < C(l +t)~^, (3.37) 

we are not able to prove convergence rate of v directly. This is because now v satisfies a 
Navier-Stokes equation which is much more complicated than the heat equation for the 
temperature variable in phase-field systems. As a result, one cannot easily obtain relation 
between ||Vw|| and v t (in certain possible norm) from the equation ( 11.11) itself. However, 
it is possible to achieve our goal by using the idea in [30] where we use higher-order 
energy estimates and construct proper differential inequalities (cf. also [5,29]). Besides, 
the convergence rate of d in higher order norm can be proved simultaneously. 

The steady state corresponding to problem (|l.ll) - (ll.5p satisfies the following system: 

^oo • V^oo - uA Voo + VPoo = -V- [VdooO Vd oo + (Arf 0O -/(d oo ))(8)d oo ], (3.38) 



V-^oo = 0, (3.39) 
^oo • V4o - doo ■ V^oo = Adoo-fidoa), (3.40) 
v oa (x+ ej) = Voo (x), doo(x + d) = d^x). (3-41) 

Lemma 13721 implies that all limiting points of system ( ll.ll) - (11.5l) have the form (0, doo) G S. 
As a result, system fl3.38j) (I3.41f) can be reduced to 

VPoo + V 0^^-) = -Vdoo-Ad^, (3.42) 
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-A^ + /(<£«,) = 0, (3.43) 
with periodic boundary condition for d^. In (I3.42p . we use the fact that 

V ■ (Vd^ VcU = V (^^-] + Woo • Adoo. 



Subtracting the stationary problem (13.42h f T3.43h from the evolution problem fll . II) f TT75ll . 
we obtain that 

/|W| 2 \ /|Woo| 



v t + v- Vv-vAv + V(P-P 00/ . 

= -X7 -[(Ad- f(d))®d}-Vd- Ad + Vdcc- Adco, (3.44) 

V-v = 0, (3.45) 

d t + v -Vd-d-Vv = A(d-d oo )-f(d) + f(d 0O ). (3.46) 

Multiplying ( I3T441) by u and ( I3T46D by -Ac/ + f(d) = -A(d - d^) + f(d) - /(<*«,), 
respectively, integrating over Q, and adding the results together, we have 

d A„ .., 

-T- - W M+ -II VCt - Vflnnir + / l-MCX) - -Mflno) - r(CXno)(CX - (Xno)ICX^ 



~||Vd- Vc/ oc || 2 + / [F(d)-F(d 00 )-f(d 00 )(d-d 00 )]d3 
1 Jq 



dt V2 

+z ,||Vt;|| 2 + ||Ad-/(d)|| 2 
= (v, Vdoo ■ Adoo) 

= (u, Vdoo • (Ac/oo - /(doc))) + (V ■ Vc/oo, /(doo)) 

= 0. (3.47) 

Multiplying (13.461) by d — d^ and integrating in f2, we have 
1 d 



, ) ^-||d-d 00 || 2 + ||V(d-d 



2 -UII 2 



= -(w Vd,d-d QO ) + (d- Vw,d-d 00 )-(/(d)-/(d 00 ),d-d 00 ) :=/!. (3.48) 
The right hand side can be estimated as follows 



|2 
oo I 



< I^II^IIVdll^lld-dooll + IIV^IIIIdll^lld-doell + linOIUoclld-d 

< C\\Vv\\\\d-d 00 \\+C(\\V(d-d 00 )\\^\\d-d 00 \\^ + ||d-doo||) 2 

< e 1 ||Vt;|| 2 + i||V(d-d 0O )|| 2 + C||d-d 0O || 2 . (3.49) 

Multiplying (13.481) by a > and adding the resultant to (13.47ft . using (I3.49h . we get 

j t QlMI 2 + |l|Vd - Vdoof + - \\d- doof + j '(F(d) - F(doo))dx 

- J /(doo)(d - doo)d^ + (!/ - aeO || Vi;|| 2 + || Ad - /(d)|| 2 + ||| V(d - d^f 

< Calld-doof. (3.50) 
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On the other hand, by the Taylor's expansion, we have 



F{d) = F^) + /(doo)(d - rfoo) + f'(0(d ~ d. 



2 

oo j i 



(3.51) 



where £ = ad + (1 — a)^ with a G [0, 1]. 
Then we deduce that 



/ [F(d)-F(d 00 )dx-f(d 00 )(d-d 00 )]dx = [ fii^d-d^fdx 
Jo Jo 



< wfim^wd-d^fK^wd-d^w 2 . 



(3.52) 



Let us define now, for t > 0, 



y(t) 



i 



a 



\\Vd(t)-Vd 00 \\ 2 + -\\d(t)-d 00 \ 



2" v 711 2 
- / f(d 0O )(d(t) - d^dx 



(F(d(t))dx - F(d 00 ))dx 



Q 



v 

2a" 



In (I3.50p and (I3.53p . we choose 

a > 1 + 2C 2 > 0, £i 
As a result, 

y(f) + C 2 |M - doof > ^(|M| 2 + |M - c/oo||h0- 
Furthermore, we infer from (13.551 ) that for certain constants C3, C4 > 0, 



^(t) + c 33/ (t) < c 4 \\d - dooir < c(i + t)-T 



29 
-29 



By Gronwall's inequality, we have (cf. [29,30]) 

28 

y (t) <c(i + ty—°, Vt>0, 

which together with (I3.55P implies that 

Wv^W + Wd^-d^Wm <C(l+t)-^s, Wt>0. 



(3.53) 



(3.54) 



(3.55) 



(3.56) 



(3.57) 



(3.58) 



Step 3. In the last step, we shall prove the convergence rate in the same space where 
the initial data stay. In Section 3.1, it has been proved that, once we could obtain the 
uniform bound of d in H 2 , we are able to obtain strong convergence of d in H 2 instead of 
weak convergence. By reinvestigating the higher-order energy estimate for the subtracted 
system (13 .441) — (13.461) (cf. also Lemma l3~ll . we can obtain a further result, which provides 
the same rate estimate of (v, d) in H 1 x H 2 as (13.581) . 

Taking the time derivative of A(t), we obtain (cf. also [26, (43)]) 



1 d 
2dt 



A(t) + (u\\Av\\ 2 + \\V(Ad - f(d)\\'' 
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= (Av, v-Vv) + (Av, V • (Vd Vd)) + (Aw, V • ((Ad - /(d)) <g> d)) 
+ (V(Ad - /(d)), V(v • Vd)) - (V(Ad - /(d)), V(d ■ Vv)) 
+(Ad-f(d),f(d)d t ) 

= I 2 + ... + I 7 . (3.59) 

Noticing that we have got uniform bounds for ||u||hi and ||d||#2 before (see Lemma I3T21) . 
in what follows we estimate Ii [i = 2, 7) term by term. 

\I 2 \ = \(Av,v -Vv)\ <\\Av\\\\v\\ L oo\\Vv\\ 

< C||Au||f ||Vu||3 < C||Au||3||Vi;||3 

< e\\Av\\ 2 + C £ \\Vv\\ 2 . (3.60) 

Before estimating 7 3 , we first look at the following estimate 

||VAd|| < ||V(Ad-/(d))|| + ||/'(d)Vd|| < ||V(Ad-/(d))|| + ||/'(d)|| L 4||Vd|| L 4 

< ||V(Ad-/(d))||+C||Ad||5||Vd||3 < ||V(Ad-/(d))|| + C. (3.61) 

Then we have 

|/ 3 | = |(Aw,V- (Vd© Vd))| = |(At;,V-(Vd0Vd))-(At;,VP oo )| 
= |(Aw, Vd- Ad) - (Aw, Vdoo • Adoo)| 

< | (Av, Vd • Ad) - (Av, Vd • Ad^) | + | (Av, Vd • Adoo) - (Av, Vdoo • Adoo) | 
=: ha + hb- (3.62) 

Furthermore, 

h a = \(Av,Vd- Ad) - (Av,Vd- Ad^l 

< ||Aw||||Vd|| LO o||Ad - Adooll 

< C||Aw||||Ad- AdoollllVdll^llVAdll^ 

< C||At;||(||Ad-/(d)|| + ||/(d)-/(d 00 )||)(||V(Ad-/(d))||+C)^ 

< e||Aw|| 2 + C E (\\Ad - f(d)\\ + ||/(d) - /(d 00 )||) 2 (||V(Ad - /(d))|| + O 

< £||Ai;|| 2 + £||V(Ad-/(d))|| 2 

+C e [l + (||Ad - f(d)\\ + ||/(d) - /(doo)||) 2 ](||Ad - f(d)\\ + ||/(d) - /(d^H) 2 

< e || At;|| 2 + e||V(Ad - /(d))|| 2 + C £ (\\Ad - f(d)\\ + ||/(d) - /(d^U) 2 

< £||Aw|| 2 + e||V(Ad-/(d))|| 2 + C e ||Ad-/(d)|| 2 + C e ||d-d 0O || 2 , (3.63) 



I 3b = \(Av,Vd-Ad 00 )-(Av,Vd 00 -Ad 00 )\ 

< ||Aw||||V(d-d 00 )|| L oc||Ad 00 || 

< ||At;||||V(d-d 0O )||5||VA(d-d 0O )||5 
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< e|| Av\\ 2 + C e ||V(d - d^W ||VA(d - d^W 

< e\\Av\\ 2 + C e ||V(d - doo)||||V(Ad - f(d))\\ + C £ ||V(d - doo)||||V(/(d) - f(d 

< eWAvf + eWViAd-fmf + CcWVid-d^W 2 

+ \\(f(d)-f'(d 00 ))Vd 00 \\) 

liv^-rfoo)!! 2 

-d 00 )\\ + \\f(d)-f(d 00 )\\ L 4vd 00 \\ Li ) 
liv^-^f 



+C £ ||V(d- 

< e\\Av\\ 2 + e 
+C £ \\V(d- 

< e\\Av\\ 2 + e 
+C £ \\V(d- 

< e\\Av\\ 2 + e 

< e\\Av\\ 2 + e 



doo)||(||/'(d)V(d-doo 
\V(Ad-f(d))\\ 2 + C £ 

doo)M\\f(d)\\L-\Wd 
\V(Ad-f(d))\\ 2 + C e 

doo)lllir(OIU-||d-d 
|V(Ad-/(d))|| 2 + C ( 

|V(Ad-/(d))|| 2 + C ( 



oo||L 4 



|V(d - d^W 2 + C e ||V(d - doo)H 5 ||d - d c 

Wd-dooW + CeWd-dooW 2 . 



(3.64) 



Next, 
h + h 



= (Av, V • ((Ad - f(d)) ® d)) - (V(Ad - f(d)), V(d • Vu)) 
= -(d • VAv, Ad - /(d)) + (Ad - f(d), A(d ■ Vu)) 

= -(d-VAu,Ad-/(d)) + [ (Adi - f i )A(d j V j v i )dx 

Jq 

= -(d- VAu, Ad - /(d)) 

+ / (A^ - ^(AdjVjVi + 2V k d j V k V j v i + djAVjVi) dx 
Jq 

= -(d ■ X7Av, Ad - /(d)) + / (Adi - fi^djVjAvi) dx + (Ad - f, Ad ■ Vv) 

Jq 

+2 / (Ad i -/ i )(V fc d j V fc V J -u i )da; 

./Q 

= (Ad - /(d), Ad • Vu) + 2 / (Ad, - /^(VfcdjVfcVjUi) dx 

i<2 



(Ad - /(d), Ad • Vv) + 2 / (Ad, - /,) [(Vd,) fe V fc (V^)J dx 
(Ad - /(d), Ad • Vu) + 2(Ad - /(d), (Vd • V) • Vu) 

£ + 4 



(3.65) 



N < ||Ad-/(d)|| L4 ||Ad||||Vt;|| L 4 

< C(||V(Ad-/(d))||^||Ad-/(d)||^ + ||Ad-/(d)||)(||A^||^||V^||^ + ||V^||) 

< e(||V(Ad - /(d))f + || A*;|| 2 ) + C £ (\\Ad - f(d)\\ 2 + \\Vv\\ 2 ). (3.66) 



|/ 6 | < C\\Ad- f(d)\\ L 4Vd\\ L 4D 2 v\\ 

< C(||V(Ad-/(d))||^||Ad-/(d)||^ + ||Ad-/(d)||)||A^|| 

< s(|| V(Ad - /(d))|| 2 + || A^f) + C e || Ad - /(d)|| 2 . 



(3.67) 
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As a result, 



h + h\ = \h + h\<\h\ + \h\ 

< 2e(\\V(Ad-f(d))\\ 2 + \\Av\\ 2 ) + C £ (\\Ad-f(d)\\ 2 +\\Vv\\ 2 ). (3.68) 



h\ = |(V(Ad-/(d)),V(vVd))| 

< ||V(Ad-/(d))||(||VT;||||Vd||Lco + ||«||L-||D 2 d||L>) 

< 1 1| V(Ad - /(d)) f + C e (|| VAd|| || Vd|| || V^|| 2 + C e || At,|| || V || 

< |||V(Ad - /(d))|| 2 + C e || V(Ad - /(d))|| || V^|| 2 + C £ ||Vt;|| 2 
+£||At;|| 2 + C £ |M| 2 

< e(||V(Ad - /(d))|| 2 + ||At;|| 2 ) + C £ ||Vt;|| 2 + C e \\v\\ 2 . (3.69) 



= (Ad- /(d), /W) 

= -(Ad - /(d), f'(d)v • Vd) + (Ad - /(d), /'(d)d • Vv) 

+(Ad-/(d),/'(d)(Ad-/(d))) 
=: I 7a + I 7b + I 7c . (3.70) 



|/ 7a | < ||Ad-/(d)||||/ , (d)|| L =o|| U || L 4||Vd|| L 4 

< C||Ad-/(d)|| 2 + C||Vt;|| 2 , 

|/ 76 | < ||Ad-/(d)||||/'(d)IUoc||d|| L oc||V^|| 

< C||Ad-/(d)|| 2 + C||VH| 2 , 

|/ 7c | < ||Ad-/(d)|| 2 ||/'(d)||^<C||Ad-/(d)|| 2 . 

Therefore, 

|/ 7 | < |/ 7 a| + |/76| + |/ 7c | <C||Ad-/(d)|| 2 + C||Vt;|| 2 . (3.71) 
Summing up, we obtain that 

\j t Mt) + iy - ^)\\Avf + (1 - 5e)||V(Ad - /(d)|| 2 
< C 5 (||V*;|| 2 + ||Ad - /(d)|| 2 ) + C,{\\vf + ||d - djfc), (3.72) 

where C^,,Cq are positive constants depending on e, ||Vi> ||, ||d||#2, ||doo||.ff2. Taking 

6 E (°'^ min 't 1 ' I/ }) ' 

we have 

\jt m + i 11 ^ 1 ' 2 + l llV{Ad " mr ~ Q>m + +¥ ~ doofHl) - (3 - 73) 
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Recalling (13.551) . multiplying (13.730 by a small positive constant a\ and adding the resul- 
tant to (13.500 . we obtain that 

j t (y(t) + %A{tj) + ^fWAvf + f ||V(Ad - f(d)f 
+ {u- ae x - C 5 a x ) \\Vv\\ 2 + (1 - C 5 ai) ||Ad - f{d)\\ 2 
+|||V(d-d 0O )H 2 

< Ca||d-d 00 || 2 + C 6 a 1 (||t;|| 2 + Hd-d^H^). (3.74) 
Constants a,e± have been chosen as in (13. 54ft . Then we can take 

ai G (0, min { — - , ^— 1 J . (3.75) 



v k 2C 5 '4C 5j , 
Using the estimates (j3.57p . (13.58)) . we infer from (I3.74p that 

±(y(t) + ^A(t))+Cr{y(t) + ^A(t)) < C 7 y{t) + C 8 (IM| 2 + \\d - dj? m ) 

< C 9 (l + t)'T=2S , Vt>0. (3.76) 

Again using Gronwall's inequality, we have 

y( f ) + -2A(t)<C 10 (l + t)-^, Vt>0, (3.77) 
which together with (13.570 implies that (cf. (13.55)) ) 

./ s 2 . . 2Cin . , 28 2Co n , , ,,o 2Cin . , 2fl 

0<A(t) < y(t) + — -{l + ty—o < — -\\d-d oa \\ 2 + — -(l + ty—e 

< Cn(l + i) _T ^, Vt>0. (3.78) 

The above estimate yields 

||Vv(t)|| + ||Ad(t)-/(d(t))||<C(l + t)-^, Vt>0. (3.79) 

Recalling (l3~25l) . it follows from (l3~79l) that 

||Ad(t) - Ad^H < C(l + t)~^, Vt>0. (3.80) 

Summing up, we can deduce the required estimate (11.81 ) from (13.581 ). ( 13.790 and (13.801 ). 
The proof of Theorem 11.11 is complete. □ 

4 Results in Three Dimensional Case 

In this section we prove the corresponding results in 3-D case, namely, Theorem 11.21 and 
Theorem 11.31 When the space dimension is three, we can still set A = 7 = 1 for the sake 
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of simplicity (cf. [18,26]). However, to prove Theorem 11.21 viscosity v plays an essential 
role, which can not be neglected. The largeness of v is needed to guarantee the existence 
of the global solution. 

The following property is useful to understand the asymptotic behavior of the solutions 
to problem (0]H jD)j>. 

Theorem 4.1. For any R > 0, whenever 

||Vt;|| 2 (0) + ||Ad-/(d)|| 2 (0) < R, 

there exists a constant Eo G (0, 1), depending on v , f , Q and R, such that either 

(1) problem (I LIB - (II .51) has a unique global classical solution (v,d) with uniform esti- 
mate 

\\v(t)\\ HHQ) + \\d(t)\\ HHQ) <C, Vt>0, (4.1) 

or 

(2) there is a T* e (0, oo) such that 

£{T*) < £(0) - s , 

where 

£(t)= 1 -\\v(t)\\ 2 + 1 -\\Vd(t)f+ f F(d(t))dx. 

J Q 

Proof. Suppose that (v,d) is a weak solution of problem fjl.ll) — fll.5j) . First, we can see 
from the basic energy law that 

£(t)<£{0), Vt>0, (4.2) 

which implies 

\\v(t)\\ + \\d(t)\\ H i<C, Vt>0. (4.3) 
By a direct calculation, we have (cf. (|3.59l) . (13. 65ft ) 

^ t A(t) + (u\\Av\\ 2 + \\V(Ad-f(d)r) 
= (Av, v-Vv) + (Av, V ■ (Vd Vd)) + (V(Ad - f(d)), V(v ■ Vd)) 
+(Ad - /(d), Ad ■ Vv) + 2(Ad - /(d), (Vd ■ V) ■ Vv) 
+(Ad-f(d),f(d)d t ). (4.4) 

We remark that the above (formal) calculation is valid for the classical solutions, but it 
can be justified by a proper approximating procedure (cf. [18,26]). We now estimate the 
right-hand side term by term. Since we only know the uniform estimate of \\v || and ||d||#i, 
the estimates we get are different from those in the previous section. 

\(Av,v-Vv)\ < \\Av\\\\v\\ L oo\\Vv\\ <C\\Av\\(\\Av\\* + l)\\Vv\\ 
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Similar to (13.611) . in the 3-D case, 

||VAd|| < \\V(Ad-f(d))\\ + \\f(d)\\ L s\\Vd\\ L e< ||V(Ad-/(d))||+C(||Ad|| + l) 

< \\V(Ad-f(d))\\+C\\Ad-f(d)\\+C. (4.6) 

As a result, 

\(Av, V- (VdQ Vd))\ 
= {Av, AdVd) < ||Ad||||Au||||Vd|| L ~ 

< C||Ad||||Ai;||(VAd||f + 1) 

< C|Ad||||A v ||(||V(Ad-/(d))||f +C\\Ad-f(d)f* +C) 

< i||V(Ad - f(d))\\ 2 + ^\\Av\\ 2 + C(\\Ad - f(d)\\ 8 + 1). (4.7) 

|(V(Ad-/(d)), V(«- Vd))| 

< ||V(Ad-/(d))||(||Vi;||||Vd|U=o + |M| L oc||Z) 2 d||) 

< C||V(Ad - /(d))||[||Vi;||(VAd||l + 1) + C(||Ai;||i + l)(||Ad|| + 1)] 

< l\\V(Ad - f(d))f + pA^| 2 + C(\\Vvf + \\Ad- f(d)\\ 8 + 1). (4.8) 

\(Ad-f(d),Ad- Vv)\ 

< \\Ad- f(d)\\ L3 \\Ad\\\\Vv\\ L 6 

< (7(||V(Ad-/(d))||'||Ad-/(d)||* + ||Ad-/(d)||) 
x(\\Ad-f(d)\\ + l)(\\Av\\ + \\Vv\\) 

< i||V(Ad - f(d))\\ 2 + ^||At;|| 2 + C(||Ad - f(d)f + \\Vv\\ 2 + 1). (4.9) 
2\(Ad- f(d),(Vd- V) • Vi>)| 

< ClIAd-JlUsllVdlUe^ll^ 

< C(||V(Ad - /(d))||*||Ad - /(c# + \\Ad - f(d)\\)(\\Ad\\ + l)(||At;|| + 1) 

< i|| V(Arf - /(d))|| 2 + ^||At;|| 2 + C||Ad - /(d)|| 6 + C. (4.10) 

As before, the last term on the right-hand side can be expressed into three terms 

(Ad-f(d),f'(d)d t ) = -(Ad-f(d),f'(d)vVd) + (Ad-f(d),f'(d)d-Vv) 

+(Ad-f(d),f(d)(Ad-f(d))). (4.11) 

Then we have 

\(Ad-f(d),f'(d)vVd)\ < \\Ad-f(d)\\\\f(d)\\ L3 \\Vd\\ L e\\v\\ L ~ 
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< C\\Ad-f(d)\\(\\Ad\\ + l)(\\Avp + l) 

< ^\\Av\\ 2 + C(\\Ad-f(d)\\f + 1), (4.12) 

\(Ad-f(d),f'(d)d-Vv)\ < ||Ad-/(d)||||/'(d)IU3|MIU-l|V^|U6 

< C7||Ad-/(d)||(||Ad||+l)(||At;|| + l) 

< ^\\Av\\ 2 + C(\\Ad-f(d)r + l), (4.13) 

\(Ad-f(d)J'(d)(Ad-f(d)))\ 
= ||/'(d)|U3||Ad- f(d) \\ls 

< C(\\V(Ad-f(d))\^Ad-f(d)f* + \\Ad-f(d)\\Y 

< ±\\\7(Ad-f(d))\\ 2 + C\\Ad-f(d)\\ 2 . (4.14) 

Summing up, we can conclude that 

^(t)<a(A(t) 4 + l), (4.15) 

where C* is a constant that only depends on u, /, Q, \\vq\\ and Hdollii" 1 - 
If the initial data satisfy 

A(0) = \\Wv \\ 2 + \\Ad -f(d )\\ 2 <R, 

we consider the following initial value problem for a nonlinear ordinary differential equa- 
tion: 

±Y(t)=C^Y(t)' + l), Y(0) = A(0)<R. 
dt 

We denote by I = [0,T max ) the maximal existence interval of Y(t) such that 

lim Y(t) = oo. 

L rJ - max 

On the other hand, it is easy to see that for any t G I, < A(t) < Y(t). Conse- 
quently, A(t) exists on /. Moreover, T max is determined by Y(0) and C* such that 
Tmax = T max (Y(0), C*) is increasing when Y(0) > is decreasing. We can take t = 
^T max (R,C*) > 0. Then it follows that Y{t) as well as A(t) is uniformly bounded on 
[0, to]- This fact together with the argument in [26] and Lemma [2721 implies the local 
existence of a unique (classical) solution of problem (ll.ip - (ll.5l) at least on [0, to]. 
If (2) is not true, we have 

£ (t) > E (0) - e , Vt>0. 
From the basic energy law, we infer that 

/»oo P 

/ / {v\Vv(t)\ 2 + \Ad(t)- f(d(t))\ 2 ) dxdt<e , Vt>0. 
Jo Jo 
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Hence, there exists a G [y, to] such that 

is\\Vv(Q\\ 2 + \\Ad(Q - f(d(Q)\\ 2 < ^. 

to 

Choosing e such that 

2 £ »< fl . 



min{l, u} to 

we have < R. Taking t* as the initial time, we infer from the above argument that 

A(t) is uniformly bounded at least on [0, ^j-] C [0, t* + to]. Moreover, the bound only 
depends on R, C* but not on the length of existence interval. We can extend the local 
unique classical solution step by step to infinity such that 

A(t) <C, V t > 0, (4.16) 

where C is uniform in time. The proof is complete. □ 



Remark 4.1. Theorem implies that if the energy £ does not "drop " too fast, problem 
( ll.ll) - (11.5l) admits a global unique classical solution. This assumption can be verified for 
certain special cases, which are stated in the following corollaries. 

Corollary 4.1. Let d* G H 2 (Q) be an absolute minimizer of the functional 

E{d) = ^||Vd|| 2 + / F{d)dx 

in the sense that E(d*) < E(d) for all d G H^{Q). There exists a constant a G (0, 1] that 
may depend on v , f , Q and d* such that for initial data (vo, do) G V x H 2 {Q) satisfying 

\\vo\\m + \\do — d*\\ H 2 < a, 
problem (11.11) - (11.51) admits a unique global classical solution. 
Proof. Without loss of generality, we shall assume that a < 1. From assumption 

IKIU 1 + I Mo - d*\\ H 2 < a < 1, 

we infer that 

u\\v \\ 2 H i + \\Ad ~ f(d )\\ 2 

< u\\vo\\ 2 H1 + 2\\Ad - Ad*f + 2\\f(d ) - f(d*)\\ 2 

< i^idl^ollj/i + ||^o — 

< K x . (4.17) 
In addition, since d* is the absolute minimizer of E(d), we have 

6(0) -6(t) < 6(0) - E{d{t)) < 6 (0) - E{d*) 
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< ^lhir + ^(l|Vd || 2 -||Vcy 2 )+ / F(d ) - F(d*)dx 

< -a 2 + Ca 
~ 2 

< K 2 a. 

Here K\ and K 2 are positive constants that only depend on d*, u, f (not on a). 
Take R = Ki, e = K 2 a respectively. We can apply Theorem 14.11 by choosing 

a = min j 1, ^-T max (K 1} C») min{l, z/}j . (4.18) 

The proof is complete. □ 

Corollary 14.11 implies that if the initial velocity v is small in H l and initial director 
do is properly close to the absolute minimizer d* of functional E(d) in H 2 , problem ( 11. 1R - 
(II .51) admits a unique global classical solution. However, from the proof of Theorem 14.11 
we can somewhat relax the "smallness" requirement from H l x H 2 to L 2 x H 1 . 

Corollary 4.2. Let d* G H 2 (Q) be an absolute minimizer of the functional 

E{d) = ^||W|| 2 + / F{d)dx 
2 Jq 

in the sense that E(d*) < E(d) for all d G H^{Q). For any initial data (vo,d ) G 
V x H 2 {Q), there exists a constant a G (0,1]. which depends on v, f, Q, \\vq\\h^ and 
\\d \\ H 2 such that if 

IKH + IK - d*\\ H i < a, 
problem fll.ip ~ (jl.5p admits a unique global classical solution. 

Proof. Without loss of generality, we assume that a < 1. Set 

Kl := u\\Vv \\ 2 + \\Ado - f(d )\\ 2 < oo. (4.19) 

Moreover, we have 

5(0) - £{t) < K 2 a, 

where K 2 is a positive constant that only depend on <i*, v, f (not on a). As in the proof 
of Corollary 14.11 we take R — K\ y Sq — K 2 o and choose 

a = min|l, ^-T max (K h C*) min{l, is}\ . (4.20) 

The conclusion follows from Theorem 14. 1[ Here, we note that now a depends on the norm 
of || foil// 1 and Hrfolli/ 2 while in Corollary 14.11 a only depends on d*. □ 
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In what follows, we proceed to prove the conclusions in Theorem 1.2 and Theorem 
1.3. First, we have the following result for both cases. 

Proposition 4.1. In three space dimension case, for the unique classical solution (v,d) 
obtained in Proposition ^. 1[ Corollary \4.1\ and Corollary \4-S\ it holds 

lim A(t) = 0. (4.21) 

Proof. (1) For the large viscosity case (cf. Proposition 12.11) . after refining the argument 
in [26], we indeed have the following differential inequality (the detailed calculation is left 
to interested readers): 

Lemma 4.1. We consider 3-D case. Set A(t) = A{t) + 1. For arbitrary z/ > 0, if 
v > > 0, then the following inequality holds for the classical solution (v, d) to problem 

jA(t) <-{u- M^A{t)) \\Avf- (l - ^rMt)) \\V(Ad-f(d))\\ 2 +M 2 A(t), (4.22) 

where M l5 M 2 are constants depending on f, \Q\, \\vq\\, Hdollif 1 ? M 2 may also depend on z/ . 

Based on Lemma 14.11 we can show the uniform estimate for A{t) by using the idea 
in [18]. It follows from the basic energy law that 

/t+l _ rt+l 
A{r)dr < J A(r)dT + 1 < M, V t > 0, (4.23) 

where M > is a constant depending only on ||fo||, \\d \\m- Take v large enough satisfying 

> Mi(i(0) + M 2 M + 4M) + 1. (4.24) 
Then by (I4.22p . there must be some T > such that 

i ~ M x A(t) 
v - M!^A(t) > 0, 1 L T i2 > 0, 

for all t G [0, To]. Moreover, on [0, To], 

— A(t) < M 2 A(t). (4.25) 

(JjL 

Denote T* = supT . First we show that T* > 1 by a contradiction argument. 
If T* < 1, then 

A(T m ) < A(0) + M 2 I A(t)dt < i(0) + M 2 M. 
Jo 

On the other hand, from the definition of T*, we have 

v < m&x{M*A 2 {T*),M?A\T*)}, 
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which contradicts (14.241) . 

Next, if T* < +00, (I4.23P implies that there is a t\ G [T* — |,T*] such that 

I(ti) < 4M. (4.26) 

As a result, 

A(T.) < 4M + M 2 / A(t)dt < AM + M 2 M. (4.27) 



From the definition of T*, we have 

< maxlM^i 2 ^),!/^ 4 ^)}, 



which together with (j4.27l) also yields a contradiction with (I4.24p . Hence, we have the 
uniform estimate 

1 

A (t)<^ Vt>0, (4.28) 

which implies that 

A(t) <C, V t > 0, (4.29) 

where C is a constant depending on /, \Q\, \\vqWh 1 , H^oIIh 2 - Thus, we infer from (14.221) 
that 

±A{t) = jA{t) < M 2 A(t) < a 

Due to the fact that A(t) e L x (0, 00) (cf. (pTSj) ). we can conclude that (jOTll holds. 

(2) Now we consider the near equilibrium case. When the assumptions in Theorem 14. II 
(or Corollary 14.11 / Corollary 14.21) are satisfied, we have (14.161) holds for all t > 0. Thus 
(I4.15P implies 

^A{t) < C*(A 4 (t) + 1) < C, Vi>0. 
at 

By the same argument as in (1), we obtain (14.211 ). 

The proof is complete. □ 

After previous preparations, we can proceed to prove the results in Theorem 1.2 and 
Theorem 1.3. 

Proof of Theorem 1.2 and Theorem 1.3. Based on Proposition 14.11 for both large 
viscosity case and near equilibrium case, one can argue exactly as in Section 3.1 to conclude 
that 

lim (\\v\\ m + \\d - doolie) = 0. (4.30) 



Notice that we now have uniform bounds for and ||d||#2. Then we are able to 

show the estimate on convergence rate (11 .3ft for both cases. To this end, we can check 
the argument for 2-D case step by step. By applying corresponding Sobolev embedding 
theorems in 3-D, we can see that all calculations in Section 3.2 are valid with minor 
modifications. Hence, the details are omitted here. We complete the proofs for Theorem 
IP and Theorem Q □ 
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5 Remark on the Flow with Non- vanishing Average Ve- 
locity 

We briefly discuss the flows with non-vanishing average velocity. Due to the periodic 
boundary condition (11.41) . by integration of (11.11) over Q, we get 



v(t)dx = 0, (5.1) 



d ( 1 
dt \\Q\ J Q 

which implies 

mv := T7TT / v ( t ) dx = T7T\ I v ° dx ' V t - °' ( 5 - 2 ) 

where \Q\ is the measure of Q. 

Our main results (Theorems 1.1-1.3) in this paper are valid for the flow with vanishing 
average velocity (see the definition of function space V), namely, m v = 0. In that case, 
we can apply the Poincare inequality to v G V such that ||i>|| < C||Vu||. This enables us 
to show that under the dissipations of system (11.1R — (11.51) . the velocity of the flows will 
tends to zero and the director d will converge to a steady state. 

When the non-vanishing average flow v is considered, as for the single Navier-Stokes 
equation (cf. [27]), we set 

v — v + m v . (5.3) 
The we transform problem (II . II) ( TT75T) into the following system for variables v and d: 

v t + v-Vv- vAv + m v ■ Vv + VP = -AV • [Vd Vd + (Ad - /(d)) ® d], (5.4) 

V-v = 0, (5.5) 
d t + v- Vd + m v - Vd-d- Vv = 7(Ad-/(d)), (5.6) 

subject to the corresponding periodic boundary conditions and initial conditions 

v(x + e*) = v(x), d(x + ej) = d(x), for x e dQ, (5.7) 
v\t=o = v Q (x) = v (x) - m v , with V • v — 0, d\ t=0 = d Q (x), for x e Q. (5.8) 

Introduce the new energy functional 

g(t) = l\\v\\ 2 + ^||Vdf + A I F{d)dx. (5.9) 



2" " 2 

It is not difficult to check that system (I5.4I )- (15.8I) still enjoys the basic energy law 

—£{t) = -v\\Vv\\ 2 - A 7 ||Ad-/(d)|| 2 , t>0. (5.10) 

By a similar argument, we can still prove the global existence of classical solution (v, d) 
to problem (I5.4l) - (15.8ft under the same assumptions as in Proposition 12.11 Theorem 14.11 
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Corollary 14.11 and Corollary 14.21 Moreover, we can prove the same higher-order energy 
inequalities like Lemma |3~T1 Lemma Hj] and (14.151) for (v,d). 

As far as the long-time behavior of the global solution is concerned, following a similar 
arguments in previous sections, we can conclude that 

lim (||S(*)|| Jff i + ||Ad(*)-/(d(*))||) = 0. (5.11) 

t— >+oo 

Recalling (15.31) . we infer from (I5.1ip that 

lim \\v(t) — m v ||#-i = 0. (5-12) 

However, in general we cannot conclude anything about the convergence of d like in 
Theorems 1.1-1.3. (15.111) implies that the 'limit' of d, which is denoted by d, will satisfy 
Ad — f(d) = with corresponding periodic boundary condition. Let us look at the 
'limiting' case such that v — v — m v and d = d. It follows from (II .5p that 

E-d = d t + v - Vd = 0. (5.13) 

Consequently, d is purely transported and it (i.e., d(x(X, t), t)) remains unchanged when 
the molecule moves through a flow field with velocity m v . However, the local rate of 
change d t may not be zero, since the convective rate of change may not vanish. Hence, 
in the Eulerian coordinates, or in Q, d(x,t) may change in time. As a result, there might 
be no steady state for the director. Obviously, this is different from the situation in 
the previous sections, where all the three rates of change are vanishing in the limiting 
case. We can look at a simple example. In the case of periodic boundary condition, let 
v = (1, 0) and d(x, 0) = d (x) for x G Q. We can see that in the Eulerian coordinates, 
d(x, t) (x E Q) is a periodic function in time such that for t > 0, d(x, t) = d(x, t + 1) with 
T = 1 being the period. 
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